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Abstract. - A magnetic model is proposed for describing the incommensurate / phase of 
spin-Peierls systems. Based on the harmonicity of the lattice distortion, its main ingredient 
is that the distortion of the lattice adjusts to the average magnetization such that the system 
is always gapful. The presence of dynamical incommensurabilities in the fluctuation spectra 
is also predicted. Recent experimental results for CuGeOs obtained by NMR, ESR and light 
scattering absorption are well understood within this model. 



Since the discovery of the inorganic compound CuGeC>3 [1], much work is devoted to 
the spin-Pcicrls transition, which is to be viewed as a magneto-elastic distortion induced by 
quantum magnetic fluctuations [2], [3]. It is usually observed in s = 1/2 antiferromagnetic 
(AF) chains with isotropic spin-spin interactions. For such uniform chains — the corresponding 
lattice structure defines the U phase of the system — the interplay between spins and phonons 
gives rise to a lattice distortion at low temperature (T <C J). Depending on the value 
of the applied magnetic field H , this distortion results in a new lattice structure which 
remains commensurate or becomes incommensurate. In low fields, the distortion corresponds 
to a lattice dimerisation — this structure defines the D phase. In large fields, the lattice 
incommensurability — which defines the / phase — increases with H [2] . Our purpose is to 
discuss for that / phase the properties of a model Hamiltonian able to describe the features 
observed experimentally. Comparisons with recent data obtained on CuGeOs are finally 
presented. 

As a result of the incommensurate lattice distortion, the magnetic Hamiltonian in the / 
phase should display similar incommensurate periodicities. For H > H c , where H c is the 
critical field of the first order transition between the D and / phases, we propose to describe 
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the exchange coupling by the following Hamiltonian 

H = jY^[l-6icos(qr i )}S i S i+1 (1) 

i 

where 5i is the "alternation" of the next-nearest exchange coupling J and q the wavevector 
characterizing the magnetic modulation. Higher anharmonicities could easily be accounted 
for by adding terms — S n cos(nqri) — with n odd integer — in the square bracket of (1). For 
q = tt, the modulation is commensurate and the D phase is recovered. The case 5\ = in 
(1) corresponds to the U phase. In the / phase, q ^ tt holds and we define q = d + tt where 
d evaluates the lattice incommensurability, which is related to the average magnetization per 
site (m) = N~ x J2i m i (N is the system size). In the XY limit of (1), for (m) ^ 0, the model 
is equivalent to tight-binding fermions with the creation (annihilation) operators c\ (ci) [6]- [8]. 
When an infinitesimal spin-lattice interaction g is present, instabilities occur in the system [2]. 
The largest one develops at q = 2kp = 2n(m) + tt, where the magnetic incommensurability 
d = 2ir(m) is created. Due to this instability, an energy gap A + + A_ opens in reciprocal 
space exactly at the two Fermi points as shown in fig. la). The states above and below these 
gaps stay separated since they evolve continuously on increasing g, i.e., no state can "jump" 
across the gap. This picture remains valid when interactions between fermions are taken into 
account — i.e., in the isotropic Heisenberg limit of (1). The corresponding chemical potential 
\i is smaller than the Zccman energy h = g^H, as the interactions create an internal field 
opposed to H. In the / phase, \x does not need to lie exactly in the middle of the energy gap 
(see fig. la)), i.e. A + ^ A_ in general. Three different gaps have to be considered: A + , A_ 
and A , which occur in the energy spectra of the spin fluctuations S±(q,E) and S\\(q,E) 
observed in directions perpendicular and parallel to the applied field, respectively (see figs. 
lb,c). 

Since the system is gapful, there exists an infrared cutoff and a conventional perturbation 
treatment [8] offers a reliable approach. The local magnetization m, = (c\ci) — 1/2 is evaluated 
as follows. In the XY limit, a chain of tight-binding fermions with hopping ti and local 
energies [ii is considered. For each bond i, an effective two-site Hamiltonian Hi and on-site 
and nearest-neighbor Green functions are defined via 

Hi= (^ + Li(E) U \ and f Gi-i G, ti+1 \ ={E _ H) -i (2) 

The self-energy Li{E) (Ri + i(E)) takes the half-chain to the left (to the right) of bond i into 
account. The conditions Li = tf_ 1 /(E — [ii-\ — Li-i) and Ri = tf/(E — — Ri+i) define 
recursively continued fraction representations. The expectation values (c\cj) are computed by 
the integral (c\cj) = (27r) _1 J r Gi-j(E)dE, where the complex contour I\ contains the real 
interval (— oo,/x]. With /j within the gaps (see fig. la)), this integral is easily evaluated. In 
the isotropic Heisenberg model, a Hartree-Fock (HF) treatment is used. The fermionic chain 
is defined with t; = — J/2 — s J(c|c,+i) and Hi = h — sJ(mi-i+mi + i) [8]. The renormalization 
s < 1 is assumed to be independent of d and H for small fields (h <C J). It is chosen as in 
ref. [8], such that in the D phase — the h, d — limit of our model — A = A + = A_. The 
self-consistent HF calculation is done here with full spatial dependence unlike the approach 
of Fujita and Machida (FM) [9] where only a constant Fock term is considered. In the XY 
limit, one obtains the results shown in fig. 2a. The local magnetization m, has an alternating 
part and a slowly varying one, which is spatially modulated with periodicity L/2 (L = 2n/d). 
With interactions between fermions — i.e. in the isotropic Heisenberg limit of (1) — these two 
important features are kept. However, the rm become (alternately) negative — see fig. 2b) - 
and their amplitude is much larger. For the same model (1), we performed also density-matrix 
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renormalisation group (DMRG) calculations on finite chains [11]. The results confirm our 
renormalised HF approach very well. Note that such antiparallel local magnetizations are 
generic to inhomogeneous AF chains [10]. 

From the fermionic dispersion sketched in fig. la), a representation of the spin fluctuation 
spectra in the / phase can be proposed. Based on arguments similar to those used for uniform 
chains [7], [12], one is led for Su(q,E), and S^(q,E) to figs, lb) and c), respectively. As 
in [12], one has to distinguish between S_ + (q,E) = (S~(q)[E - (H - E )}- 1 S+{q)} and 

(q,E) = (S + (q)[E — (H — Eoj^^-S^iq)} due to the broken spin rotation symmetry. 

Accordingly, S±(q,E) is defined as S±(q,E) = [S-+(q, E) + S + -(q, E)]/2. The description for 
the spins is obtained through = c\ exp[i7r^~]_ 00 (l — cjci)]. In mean-field treatment, this 

phase factor can be written as « exp[i7r^J~i 00 (l — (cjcj))]. With (cjcj) w 1/2 + (to) , it is 
seen to lead to a wave vector shift {-K — d)/2 [7]. Hence, for a given value of d, the lowest energy 
branch of S\\(q, E) develops an incommensurate feature at q = it — d, where the energy gap A 
occurs also (see fig. lb). In general, Ao < A + + A_. The same gap Ao appears at q — 0, while 
a Zeeman shift h develops at q — it (as in the U phase [12], it implies \i = ft/2). For S±(q,E) 
(fig. lc), similar spectra are obtained but shifted by tt, with, however, the occurrence of the 
gaps A+ and A_ at q = d. The Zeeman shift is found now at the center of the Brillouin zone 
(q = 0). In the inset of fig. lb),c), the gap A + + A_ is calculated as a function of (to) for 
the values J = 120K, Si — 0.12 (corresponding to CuGeC>3 [3]) and S3 = — 0.075i (see below). 
Renormalized HF (solid line) and DMRG (circles) results are in qualitative agreement. By 
DMRG, it was also verified that A + and A_ are indeed different, but of the same order of 
magnitude as in fig. lc). The same approaches correctly predict for the D phase [8] a triplet 
excited state giving a unique gap for H = 0: A = A + = A_ ~ 0.37 J ~ 44K [13]. For H ^ 0, 
a Zeeman splitting occurs and three branches are observed as in CuGe03 [3]. 

The lattice incommensurability in the I phase of CuGeOs has been established by X-ray 
measurements [4]. A small anharmonicity has also been observed [5]. However, the intensity 
ratio between the first and the third harmonic super-lattice peaks is small (I3/I1 ~ 1/200 [5]), 
which let us expect only small anharmonicities in the exchange coupling (63 /Si = — \fl3jl1 ~ 
—0.07). Frustrating next-nearest neighbor interactions J2 (> 0) have also been considered 
by DMRG. As shown in fig. 2b) for a = J2/J = 0.35 [15] and for the same gap in the D 
phase (E ~ 44K) the results do not change very much. For the local magnetization TOj in the 
/phase, we refer to recent high field nuclear magnetic resonance (NMR) measurements [16]. 
Above H c (~ 12.5T in CuGe0 3 ) and below T sp — T sp (~ 10K for H - 15T) is the critical 
temperature of the second order transition between the U and / phases — a distribution of the 
TO, develops in the crystal. Below 4K, the NMR spectra (fig. 4 in [16]) become practically T 
independent. In other words, the distribution of the local magnetization becomes quasi-static, 
and our T = calculation can be used. In [16], the NMR data were analyzed within the FM 
soliton model [9]. In that model, the lattice distortion refers explicitly to the sine-Gordon 
equation (which predicts too large 83/81 ratios [5]) and the spin system is treated in the XY 
limit except for a constant Fock renormalization of the coupling (for the D phase, the XY limit 
predicts only a doublet state, not a triplet as observed in CuGe03 [3]). The FM approach is 
therefore inadequate for CuGe0 3 . As shown in fig. 2b), the isotropic Heisenberg limit of (1) 
predicts a strongly negative part for TOj, which seems to contradict the NMR data [16] where 
the distribution of the local field is observed to be mainly positive and to have a much smaller 
amplitude (AS Z = S^ nlix — 5^ in ~ 0.065). This important reduction of the experimental 
value of AS Z can be explained by the zero-point motion of the so-called phasons [17]. Such 
gapless excitations are linked to the possibility for the incommensurate distortion to slide 
along the chain without energy cost. This non-adiabatic effect induces "oscillations" of the 
magnetic pattern, which reduce the alternating component of nn and result in an averaging 
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over adjacent sites as m® ff = (1 — 2j)mi + 7(m i+1 + mj_i). In the present work, 7 is used as 
a fit parameter, though a quantitative evaluation of the averaging, which is to be presented 
elsewhere, gives correct orders of magnitude. From the diamond data in fig. 2b) (a = 0, 
i.e. no frustration) one obtains for 7 = 0.20 the effective local magnetization depicted in 
fig. 2c). From the square data (a — 0.35) of fig. 2b) almost identical results (not shown) 
are obtained for 7 = 0.19. The averaging restores a mainly positive distribution, which is 
remarkably similar to that obtained in the XY limit (compare figs. 2a) and c)). Following the 
same procedure as in [16], we obtain from fig. 2c) the NMR lineshapes displayed in fig. 3. A 
reasonable agreement is achieved for 5i = 0.12, a result consistent with the value determined 
previously for CuGeC>3 [14]. For the intrinsic damping of the NMR line, we took a = 0.035T. 
This damping partly results from the experimental procedure since, in [16], the NMR signal 
was recorded as a function of H. During the field sweep, the incommensurability d is varied, 
changing continuously the distribution of the m, . This could explain the apparent "training" 
observed at low and high fields. 

Concerning the spin dynamics in the / phase, we refer to inelastic Raman scattering (IRS) 
[18], [19] and electron spin resonance (ESR) [20] measurements performed in CuGeC>3 well 
above H c . Concerning IRS, for any field values a striking result is obtained: a peak at E <~ 230 
cm -1 is observed in all the three phases [18], [19]. Such a peak occurs when the density of 
states of elementary excitations becomes large. In the D phase, it corresponds to (twice) the 
maximum of the triplet S = 1 elementary excitations (fig. la) in [18]), and in the U phase, 
to (twice) the maximum (E max ) of the low energy spinon branch (fig. lb) in [18]). The peak 
observed in the / phase (fig. 2 in [19]) is noticeably similar to the peak in the U phase. This 
strongly supports the representations proposed in figs, lb) and c) where both the U and I 
phases are represented: at high energy, the excitation branches behave similarly, with the same 
-Emax and therefore the same density of states. In an ESR measurement, the q = mode of 
S±(q, E) is directly probed. In high fields, a remarkable behavior of the ESR signal has been 
reported [20]: no change on its position and lineshape occurs when the [/-/transition line is 
crossed (fig. 1 in [20]). This result supports well our conjecture presented in fig.lc) since the 
q = mode undergoes the same Zeeman shift h in the two phases. At T sp , the only dynamical 
changes occur in the low energy part of the spectra, with the opening of the energy gaps A + 
and A_. 

In conclusion, a magnetic Hamiltonian is proposed to describe the properties of the /phase of 
a spin-Peierls system [21]. The lattice distortion is considered to result mainly in a modulation 
of the exchange alternation. Both the local magnetization distribution and the dynamical 
properties known at present for the / phase of CuGeOs are well explained by this model. It 
predicts the adiabatic / phase to be characterized at low energy by two important features: 
the presence of a dynamical incommensurability (as in the U phase) and the opening of energy 
gaps (as in the D phase). Varying the external field tends to reduce either A + or A_. The 
whole spin-lattice system, however, responds by changing its global modulation in order to 
minimize its total energy so that eventually none of the gaps closes. 
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Fig. 1. - a) Fermionic dispersion for Hamiltonian (1). The Fermi energy is denoted by /i; b) 
and c) dispersions of the spin excitations in Sn(q, E) and S±(q,E) for the / and U phases 
(solid and dashed lines, respectively). For h <C J,Em&x/J ~ (tt/2). Inset: dependence of the 
energy gap A + + A_ as a function of the incommensurability d (solid line), for Si = 0.12 [14], 
(5 3 = — 0.07<5i and a = (see text). The dots depict DMRG results for a 66 site chain. The 
dashed line is twice the gap value (0.37 J) determined for the D phase [14]. 
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Fig. 2. - a) Local magnetization versus site number in the XY limit of Hamiltonian (1) for 
L(= 2ir/d) ~ 66 (as in [16]), Si = 0.12, (5 3 = — 0.07<5i and a = 0; b) diamonds: as previously 
but in the Heisenberg limit of (1) (DMRG provides almost identical results); squares: DMRG 
results (same gap but a = 0.35, 5i = 0.033, (5 3 = — 0.07<5i for a 66 site chains ; c) the effective 
local magnetization in the isotropic Heisenberg limit for 7 = 0.20 (see text). In a), the dashed 
and dot-dashed lines correspond to the FM model [9] calculated for the value k = 0.9 (k is the 
modulus of the elliptic Jacobi function; the offset occuring in the FM description is determined 
by (to) — 1/L and L = 66. Full symbols: odd sites; empty symbols: even sites; curves through 
symbols are guides to the eye. 




Fig. 3. - NMR lineshapes evaluated from the distribution of the local magnetization displayed 
on fig. 2c), for different values of Si and S3 = — 0.07<5i. The internal NMR damping is 
a = 0.035T (see text). 



